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Abstract 

We show that certain infinitesimal operators of the Lie-point sym- 
metries of the incompressible 3D Navier-Stokes equations give rise to 
vortex solutions with different characteristics. This approach allows 
an algebraic classification of vortices and throws light on the alignment 
mechanism between the vorticity u> and the vortex stretching vector 
S co, where S is the strain matrix. The symmetry algebra associated 
with the Navier-Stokes equations turns out to be infinite-dimensional. 
New vortical structures, generalizing in some cases well-known con- 
figurations such as, for example, the Burgers and Lundgren solutions, 
are obtained and discussed in relation to the value of the dynamic 

angle <j) = arctan ^ w _^^ . A systematic treatment of the boundary 

LO -id 

conditions invariant under the symmetry group of the equations un- 
der study is also performed, and the corresponding invariant surfaces 
are recognized. 

1 Introduction 

Let us consider the incompressible 3D Navier-Stokes equations 

Ai = ut + Tt- VTt + Vp- vV 2 Tt -J* = 0, (1.1) 
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A 2 = V-u* = 0, (1.2) 

where u= (ui(x,y,z,t),U2{x,y,z,t),u 3 (x,y,z,t)) is the velocity field, p = 

p(x,y,z,t) the fluid pressure, v the viscosity coefficient, and /=/ (x,y,z,t) 
is an external force. 

Generally, in the studies of turbulence flows it is important to recognize 
the generation mechanisms and evolution of vortical structures. In this task 
both analytical and numerical methods are of help. For example, Tanaka and 
Kida [1] showed that high vorticity with relatively low strain rate corresponds 
to the vortex tube, while high vorticity with comparable strain rate to the 
vortex sheet. Wray and Hunt [2] presented an algorhitm for classification of 
turbulent structures where the flow field is classified into four regions (eddies, 
convergence, shear and stream regions). 

Interesting results for the 3D Navier-Stokes equations, adressed to this 
line of research, have been found recently by Galanti, Gibbon and Heritage 
(GGH) [3] and by Gibbon, Fokas and Doering (GFD) [4]. The first au- 
thors investigated the mechanism of the vorticity alignment in the Navier- 
Stokes isotropic turbulence. They started from the results obtained in [5], 
according to which the vorticity vector oo aligns with the intermediate eigen- 
vector of the strain matrix S. This problem was tackled via the intro- 

- A A 

duction of the variables a = ^ = £ ■ S £ (the "stretching rate") and 

Ul-U) 

_> — -» A A A - 

X— ^ b Jf = £ AS £ where £ — — . Two differential equations were derived 
for a and \ which were exploited to discuss the vorticity alignment in terms 
of the angle <j>(x, y, z, t) between u and S ou defined by tan = ^ A5 ^ = ^ 

uj-Slu 



(x — yX-X). On the other hand, GFD observed that stretched vor- 
tex solutions of the 3D Navier-Stokes equations, such as the Burgers vor- 
tices, are characterized by a uni-directional vorticity stretched by a strain 
field decoupled from them. Then, they showed that these drawbacks can 
be partially circumvented by searching classes of solutions of the type u— 
(ui(x, y, t), u 2 (x, y, t), j(x, y,t)z + W(x,y,t)). In such a way the equations 
for the third component of vorticity u 3 and W decouple. Generalizations 
of the Burgers vortex type are constructed and various solutions for W are 
discussed. 

Following the spirit of the above mentioned works, here we show as vor- 
tex solutions of the incompressible 3D Navier-Stokes equations (1. 1)— (1.2) 



2 



are generated by their symmetry infinitesimal operators. Using a group- 
theoretical approach, new classes of exact solutions are found. Some of them 
contain as special cases the Burgers vortex and the shear-layer solutions, 
which are produced by two different infinitesimal operators. By an algebraic 
point of view, equations (1. 1)— (1.2) are characterized by a symmetry alge- 
bra containing arbitrary functions. This fact indicates that this algebra is 
infinite-dimensional. Subalgebras isomorphic to the algebra of the Euclidean 
group are present. 

An important feature of the technique applied is that the symmetry alge- 
bra turns out to be not trivial only if the external force / is derivable from 
a potential function <p(x, y, z, t), i.e. 

/= Vy?. (1.3) 

The problem of the boundary conditions is also considered. This is a funda- 
mental aspect of the symmetry analysis, because if the boundary conditions 
are also invariant under the symmetry group, then an invariant solution is 
the unique solution of the system under investigation [6]. In this way, we 
have found interesting invariant surfaces which play the role of invariant 
conditions. 

The paper is organized as follows. In Section 2 the infinitesimal operators 
of the (Lie-point) symmetries admitted by Eqs. (1. 1)— (1.2) are obtained, 
while the corresponding group transformations are described in Section 3. 
Section 4 is devoted to the determination of the reduced equations coming 
from the generators of the symmetry subgroups. Some of these equations 
can be exactly solved and give rise to solutions of the vortex type which 
generalize well-known configurations as, for example, the Burgers shear- 
layer and the Burgers vortex. Section 5 deals with the invariance properties 
of the boundary conditions. Finally, in Section 6 and in Appendices A and 
B some comments and details of the calculations are reported, respectively. 

2 Group analysis 

To the aim of finding the symmetry group of Eqs. (1.1)— (1.2), let us introduce 
the vector field 
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V = ^id x + &d y + &d z + Udt+ E </>jdu s + 4 «9 P , (2.1) 

j'=i 

where £j and 0j are functions of x, y, z, t, u±, u 2 , 113, p, and d x = d/dx, and so 
on. 

We remind the reader that for a given system of differential equations 
Aj = 0, defined on a differentiable manifold M, a local group of transfor- 
mations G, acting on M, is a group of symmetries of Aj = if and only 
if 

pr {n) V[Aj] = (2.2) 

whenever Aj = 0, for every generator V of G (see [7], [8] and [9]). 
In the case of Eqs. (1.1)— (1.2), the condition (2.2) becomes 

pr (2 V[Aj] = 0, (2.3) 
where pr^V denotes the second prolongation 

pr (2) V = V + <f>*d uix + <j>\d uiy + 4>ld uiz + cf>\d uit + ^d UXxx + <t>\ y d uiyy + 

<t>?du lzz + <P x 2 d U2x + <f) y 2 d U2y + <p z 2 d U2z + 0*<9 M2t + (f>2 X d U2xx + (f>2 y d U2yy + 
<P Z 2 Z 9u 2zz + <Ptd U3x + 4>ld u . Ay + <Pld U3z + <p\d uzt + <pfd UZxx + <j>f d U3yy + 

4>l z d U3zz + 4>ld Px + 4>\d Py + <f>ld Pz , (2.4) 

with 

<f>- = D a {^- £ E (2-5) 

i=l i=l 

j = 1, 2, 3, 4, = §j^, = fff> Xi and a stand for x ' y ' z > *' and x ' ^' z ' 
i, xx, 1/1/ , zz, respectively, .D being the derivation operation. 

Equation (2.3) provides a set of differential equations, the so-called deter- 
mining system, which allows us to obtain the coefficients £j,(f)j- From these 
we infer the infinitesimal operator generating the symmetry group of Eqs. 
(1.1)-(1.2), namely 

V =E V it (2.6) 

i=i 
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where 

V 1 = gd x + gd ul + (-'gx + gh)d p , (2.7) 

V 2 = hd y + hd U2 + {-hy + hf 2 )d p , (2.8) 

V 3 = rd z + rd U3 + (-rz + rf 3 )d p , (2.9) 

V A = kd p , (2.10) 



3 

V 5 = a[xd x + yd y + zd z + 2td t - u i9 Ui 



i=i 



+ (-2p + xh + Vh + *h + 2 f fidx + 2t f htdx)d p \, (2.11) 

V 6 = b[yd x - xd y + u 2 d Ul - u x d U2 + {yj\ - xf 2 )d p ], (2.12) 
V 7 = c[zd y - yd z + u 3 d U2 - u 2 d U3 + (zf 2 - yf 3 )d p }, (2.13) 
V s = d[zd x - xd z + u 3 d Ul - uxd U3 + (zf x - xf 3 )d p ], (2.14) 

V 9 = e[d t + (f f lt dx)d p ], (2.15) 

g(t), h(t), k(t), R(t) are arbitrary functions of time of class C°°, and a, b, 
c, d, e arbitrary constants (dots mean time derivatives). The infinitesimal 
operators (2.7)-(2.15) represent the generators of the Lie-point symmetries 
of the 3D Navier-Stokes equations (1.1)— (1.2). 
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3 Group transformations 



The integration of the infinitesimal operator (2.6) enables us to find the finite 
transformations leaving the equations (1.1)— (1.2) invariant. We have that the 
linear combination W — V\ + Vi + V3 gives rise to a transformation to an 
arbitrary moving coordinate frame of the type 

: (x, t, u,p) — ► (x +e a, t, u +e a,p—e x ■ a — e 2 a ■ a + a • f de) 
a ' 2 Jo 

(3.1) 

where a = (g(t), h(t),r(t)) and e is the group parameter. 
V4 is the generator of the pressure change 

G p : (x,t,u,p) — > (x,t,u,p + £ k(t)) , (3.2) 

while V 5 is the infinitesimal operator for the scale transformations. The 
vector fields V 6 ,V 7 ,V 8 yield both the space and the velocity field rotations. 
Finally, Vg produces the time translations together with an adjustment of 
the pressure, i.e. 

G : (x,t, u,p) — > (x, t + e,u, p + ip(x,y, z,t + e) — <p(x,y, z,t)j , (3.3) 

where ip is the potential function linked to the external force by (1.3). Then, 
we can " absorbe" / in p by re-defining the pression p as 

p' = p - (p( x ,y,z,t). 

Without introducing new symbols, we shall put formally /= into Eqs. 
(1.1)— (1.2) and into the expressions (2.7)-(2.15) for the generators. 

4 Reduced equations and exact solutions 

By exploiting the generators Vj of the Lie-point transformations (see (2.7)- 
(2.15)), one can build up exact solutions of Eqs. (1.1)— (1.2) via the symmetry 
reduction approach. This allows one to lower the order of the system of 
differential equations under consideration using the invariants associated with 
a given subgroup of the symmetry group. In the following we present some 
reductions leading to exact solutions of Eqs. (1.1)— (1.2) of physical interest. 
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4.1 Case i) 

Let us take the vector field V\ defined by (2.7). A set of basis invariants Ij of 
the related subgroup can be determined from the finite group transformations 



x = x + eg(t), y' = y,z' = z,t' = t, 
u[ = ui+e9,u' 2 — u 2 ,u' 3 — u s , 

p' = p-\e'9{2x + eg). (4.1) 



We obtain 



9 r t t . 9 2 



h =V,h = z, h = t,l± = ui x, h = u 2 , h = u 3 ,I 7 = p + —x . (4.2) 

9 2g 

By means of the choice of variables 

g 

Ui(y,z,t) = h = ui - -x,U 2 (y,z,t) = I 5 = u 2 ,U 3 (y,z,t) = u 3 , 

n(y,z,t) = I 7 =p+^-x 2 , (4.3) 
the system (1.1)— (1.2) is cast into the reduced form 

U lt + -Ui + U 2 U ly + U 3 U lz - u(U lyy + U lzz ) = 0, (4.4a) 

U 2t + U 2 U 2y + U 3 U 2z + n y - u(U 2yy + U 2zz ) = 0, (4.4b) 



U 3t + U 2 U 3y + U 3 U 3z + 7i z - u(U 3yy + U 3zz ) = 0, (4.4c) 

U 2y + U 3z + - = 0. (4.4d) 
9 

A particular solution of Eqs. (4.4a)-(4.4d) can be found as follows. Let 
us put [10] 

z 2 y 2 

U 2 = hy, U 3 = a - kiz, n = + faaz - k\— + r(t), (4.5) 
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where ki,k2,<J,g are constants, and r(t) is an arbitrary function of time. 
Substitution from (4.5) into Eqs. (4.4a)-(4.4d) yields the linear equation 

U u + hyU ly + (a- hz)U lz - u(U lyv + U lzz ) = 0. (4.6) 

At this point it is convenient to look for a solution of Eq. (4.6) of the 
form 

U 1 = Y(y)T(z)$(t). (4.7) 
Then, Eq. (4.6) can be written as 

Y Y T T $ + 

hy^ -»-f = (hz - a) t± + u 1 -^ - | = G, (4.8) 

G being an arbitrary constant. On the other hand, Eq. (4.8) tells us that 

(hz - + z/% = G + $ = H, 

where H is an arbitrary constant. To summarize, the functions Y, T, $ obey 
the ordinary differential equations 

$ t = (H - G) $, (4.9a) 



vY m - k iy Y y + GY = 0, (4.9b) 

z/T zz + (Jfe^ - a)T z — HT — 0. (4.9c) 
Equations (4.9a)-(4.9c) afford, respectively, the general solutions 

$ = ciexp[(#- G)t], (4.10a) 

Y = ^"S? 2- -£> + " C3M( 2 " S? 2' (410b) 

if 1 -kAz-j-) 2 a 1 H 3 -ki(z - f~) 2 

(4.10c) 
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with ci, C5 arbitrary constants. M is the Kummer function defined by [[F] 

M a, /3, s = 1 + -2 + )-f - + ... + )-f - + 4.4 

/? (/3) 2 2! (/?) n n! 

with 

(a) n = a(a + l)(a + 2). ..(a + n - 1), (a) = 1. (4.5) 
Consequently, from (4.3) we obtain 

m r-Vrr , ,/ G 1 A^?/ 2 . „,,1 G 3 kill 2 ., 

= c ie '»- c "{[c 2 M(-^, j, -£) + S c 3 M(- - _, -, -X)] x 

if 1 -£4(2 -f) 2 a 1 H 3 -kAz-ff 

(4.11a) 

u 2 = k iy , (4.11b) 



■u 3 = cr — kiZ, (4.11c) 

2 2 

p = -kf- + fc lt 7z + r(t) - kf-. (4.11d) 

It is noteworthy that at least in some special cases, Eqs. (4.11a)-(4.11d) 
reproduce interesting solutions which can be interpreted by a physical point 
of view, such as the Burgers shear-layer and other solutions. This aspect is 
discussed below. 

4.1.1 The Burgers shear-layer and other solutions 



We observe that Eqs. (4.11a)-(4.11d) give rise to static solutions for H 
G = A. We shall distinguish two cases: a) A = 0, and b) A 7^ 0. 

Case a) 

From (4.1a) we obtain 
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«, = h + o.J—^-^,,)^^-*^,)], (4.12) 

where for simplicity <r = and we have used the properties M(0, b, <;) — 1 
and [11, p. 509]] 



erf ? = 4= [' exp(-9 2 )d8 

\/7T JO 



with 



We remark that by choosing k± = —7 (7 > 0) and C5 = in Eq.(4.12), 
we get the Burgers shear-layer solution |3|. [| 

where 

= = A^etf(^y) + B, (4.14) 
with A and B arbitrary constants. 

Now let us recall some basic quantities appearing in the study of tur- 
bulent flows. They are the vorticity 00= VA u, the strain matrix S whose 
elements are defined by S%j = \{u%j + where Uij = |^ and u^i = 
the energy dissipation SV,-SV/ and the enstrophy UiUJi (the summation under 
repeated indices is understood) [11]— [12]. Furthermore, in the study of the 
vortex alignments, it is important the dynamic angle <p(x, y, z, t) defined in 
the Introduction. 

The vorticity corresponding to the solution (4.14) is uo= (0,0,cj 3 ) t with 

lu 3 = -Y\y) = -A^expi-^f), (4.15) 
y 7 Iv 



while the strain matrix S is given by 
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S= \Y> 
\ 



/ 



-7 
7 / 



(4.16) 



The energy dissipation reads 



SySij = 2 7 2 + V 2 , 



(4.17) 



which tells us that the dissipation is high when the vorticiy is so [4]. We 
observe that the vorticity vanishes for large values of y 2 . We notice also that 
the pression p satisfies the equation 



which is just the equation (66) shown in [3]. Within our framework, p arises 
as a consequence of the reduction procedure involving the generator V\. The 
dynamic angle related to the vortex solution (4.13) turns out to be zero, 
so that the vectors oo and S to are parallel. 

From formula (4.12) we deduce that for fa > 0, by switching off the 
constant C3 a shear-layer of the Burgers type of the form u= (Z(z), —jy, 7^) T 
emerges, where now 



Therefore, (4.12) can be interpreted as a more general vortex configuration of 
the Burgers shear-layer type in which both the variables y and z are involved. 

Case b) 

Formula (4.11a) can be exploited to find other interesting static solutions 
which can be considered as generalizations of the Burgers shear-layer. An 
explicit example corresponds to the choice A = —k\ and a = c 3 = c 5 = 0. 
We have 



dp 



7T = ~^ z > 
oz 




(4.18) 




(4.19a) 



u 2 = fay, 



(4.19b) 
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u 3 = -kiz, (4.19c) 

A being a constant. 

The vorticity oo has cu 1 = and the remaining components different from 
zero. Precisely 

u 2 = -^zeM^(y 2 -z% (4.20a) 

IS SjlS 

^3 = -^/exp[^ 2 -, 2 )]. (4.20b) 

\S ZiJJ 

A possible physical interpretation of this situation is that for k\ < 0, 
it represents a vortex structure where oo 2 behaves as a shear-layer of the 
Burgers type along the axis y for any fixed value of z, while in opposition to 
ou 2 , c^3 vanishes at y — 0. The evaluation of the dynamic angle for the vortex 
(4.19a)-(4.19c) leads to the formula 

tan = ? Xy - = sinh 29, (4.21) 
y 2 — z l 

where we have put y = cosh#, z = sinh#. For example, we have = for 
9 = 0, and = ±f for 9 -> ±oo. 

An example of a non-static solution 

Non-static solutions can be obtained from (4.11a)-(4.11d) assuming that 
H — G 7^ 0. An interesting example is given by 

Ul = Aexp(-|0exp(-^)y^_i(^), (4.22a) 
u 2 = - iy , (4.22b) 
n 3 = iz, (4.22c) 



P=-y0/ 2 + ^ 2 ), (4.22d) 
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which is derived from (4.11a)-(4.11d) for H = 0,G — ^, k\ — —7(7 > 0), C3 = 
0, C5 = 0, a = 0, r = 0, where / denotes the modified Bessel function. 
The vorticity is to = (0, 0, to 3 ), with 

co 3 = -u ly = J Bexp(-^)exp(-^)yi[/_ i (^) ~ h(^)}. (4.23) 

The strain matrix and the dissipation corresponding to the solution (4.22a)- 
(4.22d) can be derived from (4.16) and (4.17), respectively, via the substi- 
tution of Y' by Uiy (see (4.13)). When v and 7 are fixed, for large y we 
have 

while the total strain behaves as 

SySy ~ 2 7 2 + ^-y~\ (4.24) 

The dynamic angle corresponding to the solution (4.40)- (4.42) is zero. There- 
fore, this vortex-like structure is characterized by a vorticity which is aligned 
with S uj . 

4.2 Case ii) 

The vector field 

V 2 = hd y + hd U2 - hydj 
is related to the invariants 

h 

t,x,z,u 2 - -ry,u h u 3 . 
h 

Then, by using the variables 

Mi = Ui(x,Z,t), 

h 

u 2 = U 2 (x,z,t) + -y, 



(4.25) 

(4.26) 

(4.27a) 
(4.27b) 
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u 3 = U 3 (x,z,t), (4.27c) 

P = k- |V, (4.27d) 
we are led to the reduced system 

U u + CWi* + U 3 U lz + 7i x - v(U lxx + U lzz ) = 0, (4.28a) 

C/ 2i + UiUix + C/ 2 ^ + U 3 U 2z - ^([/^ + C/ 2 ^) = 0, (4.28b) 

U 3t + U x U 3x + C/ 3 C/ 32 + - + U 3zz ) = 0, (4.28c) 

U lx + U 3z + ^- = 0. (4.28d) 



By changing x with y and f/2 with C/i, this system can be discussed in a way 
similar to that followed for Eqs. (4.4a)-(4.4d). Of course, starting from the 
field V 3 a reduced system analogous to Eqs. (4.28a)-(4.28d) is derived. 

4.3 Case iii) 

Let us deal with the vector field 

W = V 1 + V 2 + V 3 = gd x + hd y + rd z + 9 d Ul + h d U2 + r d U3 

+ (9x+hy+r z)d p . (4.29) 
The group transformations associated with W are 

x' = x + eg,y' = y + eh, z' = z + er, t' = t, 
u[ = u\ + e 9,u 2 = u 2 + e h, u' 3 = u 3 + e r, 

V = p - £ -^(2x + eg) - £ -^(2y + eh) - £ -^(2z + er) (4.30) 
From (4.30) we deduce the invariants 

I = t,I 1 = u 1 - x-, I 2 = u 2 - y\, I 3 = z -, (4.6) 
g h r 
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Now let us introduce the variables 

U x = i u U 2 = I 2 , U 3 = h,7t = I 4 , (4.32) 

depending on the (invariant) I — t only. Then, Eqs. (1.1)— (1.2) furnish the 
reduced system 

U lt + -C/i = 0, (4.33a) 

U 2t + \u 2 = 0, (4.33b) 

U 3t + -U 3 = 0, (4.33c) 
r 

Q h r 

l + !l + -=0. (4.33d) 



These equations give 



and 



g h r 



tfi = -, Ob = = (4.34) 



(7 h r 



UxU 2 U z = const. (4.35) 
Then, from (4.31) and (4.34) we infer 

Cl , 9 c 2 h c 3 s 

ui = \- x -,u 2 = — + y-, u 3 = — + z-, (4.36) 

g g h h s s 

while the condition (4.33d) ensures that the equation V- u— is satisfied. 
The solutions (4.36) imply that the vorticity uo = VA u vanishes, so that 
the motion is irrotational. Thus, we may introduce the velocity potential 
defined by u= V0. We have 

= l(x 2 - + y 2 \ + z*-) + C -l x +%+% + O . (4.37) 
2 g h s g h s 

Of course, fulfils the Laplace equation V 2 = 0. 
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4.4 Case iv) 

The generator V5 gives rise to the finite group transformations 
x = xexp(e),y = yexp(e),z —zexp(e),t = texp(2e), u ± = U\ exp(— e), 

u 2 = ui exp(— e), u 3 = w 3 exp(— e),p =pexp(— 2e). (4.38) 

A set of basis invariants is 

2 

tC tC Ou 

£ = -,n= -,6 = —, Ai(C, V, 0) = u x x, A 2 (£, 77, 0) = w 2 a;, 
y z t 

A 3 (£, 77, 0) = « 3 :r, tt(£, 77, 0) = px 2 . (4.39) 
The reduced equations read 

-2ir + 2^7r e + r]n ri + £7^ - 2uA 1 - A 2 + (2u - 0)0 A w + 26A 1 A W - A0 2 uA we + 

2z/77(l-r/ 2 )A 1 ^+7/A 1 A lr ,-77 2 A 3 A 1 , ) -4z/^A 1 ^-z/7/ 2 (l+77 2 )A lw +2z/e(l-e 2 )Ai ? + 
f A^ - e 2 A 2 A 15 - 4^0A^ - - z/£ 2 (l + £ 2 )Ai K = 0, (4.40a) 

-2z/A 2 - A X A 2 + (2z/ - #)#A 2e + 2OA 1 A 20 - A0 2 uA 2ee + 
2z/77(l-77 2 )A 2 ^+77A 1 A 2 ^-77 2 A 3 A 2 , ? -4z/^A 2 , ?e -z/7/ 2 (l+77 2 )A 2w +2z/e(l-e 2 )A 25 + 

eA 1 A 2? -e 2 A 2 A 25 -4^A 2 ^-2^A 25 ,-^ 2 (l+e 2 )A 2K -e 2 ^ = 0, (4.40b) 

-2z/A 3 - A X A 3 + {2u - 0)0 Km + 20AiA 3 «, - 40 2 uA 3ee + 

2z/7/(l-7/ 2 )A 3 ^+77A 1 A 3 ^VA3A 3 ^-4z/^A 3 , ?e -z/77 2 (l+7/ 2 )A 3w +2z/e(l-e 2 )A 35 + 
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eA 1 A3 f -e 2 A 2 A3 5 -4^A3 fe -2z/^A3 f r,-^ 2 (l+e 2 )A3K-^vr, = 0, (4.40c) 

- Ai + 26 A w + v A lv + £A 1? - v 2 A 3v - £ 2 A 25 = 0. (4.40d) 

A special solution of this huge nonlinear system of partial differential 
equations is derived in Appendix B. It reads 

Mi = ct"5, (4.41a) 

t~ ~ t~ ~ 

u 2 = -2at~i +ci(i/t~ 1 )* exp[— — (Ac-xt~* )]-c 2 (vt~ 1 ) 2 exp[— — (4c-rrt~2)] x 

erf i[]-v~5(2c- xt~^)], (4.41b) 

_ 1 _l 

1 1 cct ^ 1 1 cct ^ 1 

■u 3 = -26t _ 2+c 3 ( Z /t- 1 )2 exp[— — (4c-a;t"2)]_ C4 ( Z /t- 1 )2 exp[— — (4c-rrt~2)] x 

erf i[^i/~5(2c- (4.41c) 
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p = Cot' 1 + r^(-x - ay -bz), (4.41d) 



where 



erf(C) = -p/ exp(-t 2 )rft, (4.42) 

V 7T JO 



2 fC 
0F , 

erfi(C) = -ierf(iC) = ^ [* exp(-t 2 )dt, (4.43) 

Vvr Jo 

and a, 6, c and Cj (j = 0, 1, 2, 3, 4) are constants. 

The vorticity corresponding to the solution (4.41a)-(4.41d) is 

cJ= (0,-^A 3e ,-^A 2 + ^A 2e ) T , (4.44) 

V Ob L 

where A 2 (9) and A 3 (6) are expressed in Appendix B. We notice that in this 
case we obtain a stretching rate coinciding with the null vector, i.e. S co= 
(0, 0, 0) T , so that the dynamic angle = arctan(^) is not definite. 
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4.5 Case v) 

Let us consider the symmetry operator V%. The corresponding group trans- 
formations are 

x = x cos £ -\- y sin e, (4.45a) 



y = —re sine + y cose, (4.45b) 

z' = z, (4.45c) 

t' = t, (4.45d) 

U\ = u\ cose + u 2 sine, (4.45e) 

C/ 2 = —ui sine + -u 2 cose, (4.45f) 

U 3 = u 3 , (4.45g) 



which leave the quantities 



p = p (4.45h) 

z,t,r = (x 2 + y 2 )K (4.46) 
n = p (4.47) 



invariant. 
Putting 



x = r cose, y = r sine, (4.48a) 

Mi = Ui cose — C/2 sine, (4.48b) 

U2 = Ui sine + U 2 cose, (4.48c) 

% = I/3, (4-48d) 



where Uj = Uj(r,z,t), Eqs. (1.1)— (1.2) take the form 
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U lt + IWir - — + U lz U 3 + 7r r - u(U lrr + — - ^ + C/u,) = 0, (4.49) 



C/ 2t + CW* + — + U 2z U 3 - z/(f/ 2rr + — - ^ + C/ 2 ^) = 0, (4.50) 

U lr + — + U 3z = 0, (4.51) 
r 

C/ 3< + War + Ws* - v(U 3rr + + U 3zz ) + 7T 2 = 0. (4.52) 

To provide an example of exact solution to the reduced system (4.49)- 
(4.52), let us look for a solution such that 

U x = — , (4.53a) 
r 

U 2 = ^, (4.53b) 
r 

where a Q , f3 are constants. Then, Eq. (4.49) yields 

1Xr = Q ° r 3 ^° ' (4.53c) 

while Eq. (4.50) turns out to be identically satisfied. Furthermore, Eq. (4.51) 
tells us that U 3z = 0, i.e. U 3 = U 3 (r,t). Consequently, Eq. (4.52) becomes 

U 3t + ^^U 3r - uU 3rr = -n z . (4.54) 
r 

The compatibility condition n rz = n zr gives 

OL — IS 

d r (U 3t + ^—U 3r - uU 3rr ) = 0. (4.55) 
By integrating (4.53c) with respect to r we get 



<4±ft 

2r 2 

F(z, t) being a function of integration. Since n z = F z , we deduce that F(z, t) 
has to be of the form 
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F(z,t)=a(t)z + b(t), (4.57) 
where a, b are function of time. Hence 

0L — V 

U 3t H U 3r - vU Zrr = -a a , (4.58) 

r 

where we have supposed a = a Q = const, b = b Q = const. From (4.56) we 
obtain the expression for the pression p, which reads 

a 2 + 8 2 

p = 7l = -^gf^ + a z + b . (4.59) 
To find U 3 , we shall distinguish two cases. Precisely: 

I) a = v. 

Then, Eq. (4.58) takes the form of the heat equation with a constant 
source, namely 

U 3t ~ vUzrr = -do- (4.60) 

A solution of Eq. (4.60) is given by 

u ' = ^-^ + ^- (4 ' 61) 

where 7 Q is a constant. Thus, collecting all the information, we find that the 
Navier-Stokes equations (1.1)— (1.2) afford the solution 

vx ~~ (a an ^ 

M 2 = ~ 2 , (4.62b) 

u 3 = U 3 (see (4.48d)), and (see (4.61)) 

v 2 + B 2 

P= —. j^ + a z + b . (4.62c) 

II) a ^ v. 
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In this case Eq. (4.58) can be transformed into the equation 



Qt — V 

ty t + — ^ r - V^ rr = 

r 



via the change of variable 



By setting 



U 3 = *(r,t)-a t. 



# = M(t)N(r), 



Eq. (4.63) is splitted into the ordinary differential equations 



(4.63) 



(4.64) 



(4.65) 



M t = -5 M, 



(4.66a) 



vN r . 



a n -v 



N r + 5N = 0, 



(4.66b) 



where 5 > is a constant. These equations can be easily solved to give 

M = M Q exp(-St), (4.67a) 



N = r^c l J,J^r) + c 2 Y,J^-r)}, 



(4.67b) 



where M , ci, c 2 are constant, « = and J M , denote the Bessel functions 
of the first and the second kind, respectively. Substituting (4.67a) and (4.67b) 
into (4.64) provides 



5 5 

u 3 = U 3 = M exp(-5t)r^[ Cl J M (y-r) + c 2 Y^J -r)\ - a a t. (4.68) 

The components Ui,u 2 and the pression p are 

a Q x — f3 y , . 

Mi = , (4.69a) 
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u 2 = , (4.69b) 

a 2 + 3 2 

p=-^gf^ + a z + b , (4.69c) 



respectively. 



4.5.1 The Burgers vortex and other solutions 



It is noteworthy that the infinitesimal operator V6 leads to the Burgers 
vortex solution [14] as a special case, namely: 

Ul = -l x -yf(r) (4.70a) 

U2 = ~y + xf(r) (4.70b) 

u 3 = -fz, (4.70c) 

where 

/ W = /o l-exp(-a r ^ (4 n) 
7 and fo are constants, and 

a = (472) 

In fact, by choosing 

U 1 = -^r-U 2 = rf(r),U 3 = 1 z, (4.73) 
the reduced equations (4.49)-(4.52) are satisfied provided that 

TiY = r f 2 (r) - — r, tt z = -7 2 z, (4.74) 
where f(r) obeys the equation 

f r . r+ f^L r + -)f r + ^f = 0. (4.75) 
22 



The general solution of Eq (4.75) can be written as 

/ (r) = /o - (/ °-^ eXp( - ar2) , (4.76) 

which reproduces just (4.71) for f\ = 0. Eqs. (4.70a)-(4.70c) are derived 
from (4.49)-(4.52), with the help of (4.73). 
Eqs. (4.74) give the pression (p = ir) : 

V = -y (z 2 + ^-02 I 1 - ex P (-*r 2 )f+af \ [Ei (-ar 2 ) - Ei (-2ar 2 )' 

(4.77) 

where Ei denotes the exponential-integral function [11] 
while the vorticity is 



u? =(0,0,uj 3 f (4.78) 
with 

u 3 = 2f + rf r = 2af exp(-ar 2 ). (4.79) 

We point out that (4.79) coincides with the expression for the Burgers 
vortex [14], in which there exists a balance of dissipation and stretching [3], 
[15], [4]. 

Now by integrating Eqs. (4.70a) and (4.70c) we easily find 

r = r exp(--t] , (4.80) 



which tells us that in the limit t — > oo, r — > 0. For r — > both / (r) given by 
(4.71) and the vorticity (4.79) turn out to be finite. 

4.5.2 The Burgers-Lundgren solution 



Another interesting solution of the Navier-Stokes equations (1.1)— (1.2) of the 
Burgers vortex type, generated by the vector field V6, is obtained assuming 
that 

Ui = ~r,U 2 =r /(r, t),U 3 = >yz, (4.81) 
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where the function f(r,t) depends on the time also. 

In doing so, substitution from (4.139) into Eqs. (4.49)-(4.52) gets 

U-v U ~ (|r + ^) f r - 7 / = 0. (4.82) 
By inspection, Eq. (4.82) is satisfied by 

1 — exp —ar 2 \-, -, — rv) 

/ (r, t) = / 1 -Alzf^zlM (4.83) 

In this case, the vorticity is time dependent and takes the form 

7? = (0,0,w 3 M)) T , 

with 



exp 



w 3 (r, t) = 2/ (r, «) + r/ r (r, t) = 2 o Lzggizlfli . ( 4 . 84 ) 

1 — exp (—7 1 ) 

We remark that this quantity is the same as the expression of the vorticity 
given by Lundgren [15, formula (16)]. When t — ► oo, 003 (r,t) tends to the 
Burgers vortex (4.71). Therefore, as it was noted in [15] one can resort 
to the assumption of Townsend [16], who considered the axial strain rate 
proportional to the root-mean-square strain rate in turbulent flow. Under 
this hypotesis, the Burgers vortex leads naturally to the Kolmogorov lenght 



7/ 



v\ 1/4 



where e is the dissipation rate per unit mass. 

4.5.3 Another vortex-like solution of Eq. (4.82) 

Let us look for a solution of Eq. (4.82) of the form 

f = A(t)exp[B(t)r 2 }, (4.85) 
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where A and B are functions of time to be determined. Inserting (4.85) into 
Eq. (4.82) yields 

A A ° 



with a — ^L. 

By choosing = l,we obtain 

/(r , t ) = fsec^(l t )exp(- 1 + e ; r ;_ 7t) ). (4.86) 

In opposition to the behavior of (4.83) and (4.84), this solution tends to 
zero when t — > oo. 

4.5.4 Evaluation of the dynamic angles for vortices related to V 6 

1) The solution (4.61)-(4.62b) has a vorticity given by u— (ui, oj 2 , 0) T , with 

"i = [-f ~ -f {vt)-i exp(-— )}y, u 2 = [--?- + -j- (z/t)"§ exp(-— )]ar, 

(4.87) 

which leads to a constant dynamic angle, namely 

= arctan^. (4.88) 

The vorticity lies on the x, y— plane and is ortogonal to r= xx + yy, i.e. 
uo ■ 7= 0. The orientation of uo with respect to S uo depends on the ratio ^f; 
for p o — 0, = 0; for f3 — v, (f) — |, while in absence of viscosity (u — > 0) 
and when /? ^ we have = | . 

2) Let us consider the solution (4.68)-(4.69b). Owing to the complicated 
structure of 113 (see (4.68)), the expression of the vorticity results too cum- 
bersome. Therefore, it is convenient to examine the asymptotic behaviour of 
u 3 for r — > 00 at fixed ji. By choosing a Q = 0, we get [11, p. 364] 
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, _. . ,,_l r . 5 IT 1 . 8 It 7T 

m 3 ~ exp(-dt)r" 2[ciCos(W-r-/x- - -) + c 2 sm(W -r - /i- - -)]. (4.89) 



Now the vorticity corresponding to the solution (4.69a), (4.69b) and (4.89) 
and the dynamic angle can be easily calculated. For simplicity, we report only 
the value of 0, which is 

Thus, the discussion on this result is similar to that performed in the case 
1). We observe that for a Q = f3 Q = 0, is not defined. 

An interesting case is represented by the choice /i — | in Eq. (4.89). In 
fact, we get 

Q 2 

u s ~ exp( — ^-t 2 )(ci sin £lt — c 2 cos fit), (4.90) 

where f2 2 = 25 and r = Formally, the quantity (4.90) can be interpreted 
as a special case of the general solution of the parametric oscillator described 
by the equation 

Q +47it q +(72 + 271 + 4 7l 2 t 2 )g = 0, (4.91) 
whose general solution is 

q = exp(-7!t 2 )(ciCOS7 2 t + c 2 sin7 2 t), (4.92) 

where the amplitude plays the role of a damping function (7! > 0) which 
tends to zero for large values of time according to a law of the Gaussian 
type. 

3) The dynamic angles for the Burgers and Burgers-Lundgren vortices, 
together with the vortex solution with f(r,t) given by (4.86), are zero. 
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5 Invariance of the boundary conditions 



The study of the invariance properties of a system of differential equations 
with suitable boundary conditions plays a fundamental role in the description 
of realistic models. 

Generally, boundary value problems can be treated more simply for ordi- 
nary differential equations. Indeed, a symmetry analysis allows one to reduce 
a given boundary value problem for an ordinary differential equation to a 
boundary value problem for a (reduced) equation of lower order. In the case 
of a partial differential equation, by a group point of view, a boundary value 
problem is solved when one can determine a solution which is invariant, to- 
gether with all the boundary conditions, under the action of the infmitesiaml 
generator of the symmetry group. However, for a linear partial differential 
equation one can resort to less restrictive conditions [8, p. 215]. 

In this Section we present a systematic investigation of all the surfaces 
which can take the meaning of invariant boundary surfaces in a boundary 
value problem for the Navier-Stokes (1.1)— (1.2). 

To this aim, we recall that if V denotes a generator of symmetry trans- 
formations, a basis of invariants for V is obtained from the equation 

VI = 0. (5.1) 

9 

Now let V = J2 y% be the vector field generating the set of the Lie-point 

i=i 

transformations of Eqs. (1.1)— (1.2). Using the method of characteristics, Eq. 
(5.1) can be written in the form 



dx dy dz 

ax + by + ez + g Q ay — bx + cz + h Q az — cy — ex + r Q 

dt dui du 2 

2at + d —au\ + bu 2 + ew 3 — au 2 — bu x + cu 3 
du 3 dp 



—au 3 — cu 2 — eu\ k Q — 2ap 

Here we have put g = g Q = const, h = h Q = const, k = k Q = const. By 
solving the system (5.2), we can determine the characteristic manifolds A = 
const, x — const, ip = const in the coordinate space (x,y,z,t). To this aim, 
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first let us analyse the general case. Subsequently, some special interesting 
case will be examined. 



5.1 The general case 

The system (5.2) can be put into the form 

die 



dr 



= A-x + r (5.3) 



where: 





I a 


A = 


u 






x \ 




I, 


V 








(5.4) 



x = \ y , r = \ h (5.5) 



r= lln^ + l). (5. 6 ) 




2a V d 

Let us distinguish the following cases. 

Case I) a ^ 0. 

Eq. (5.3) admits the solution 

~~x = e Ar ~xl + xp^ (5.7) 
where ~xo is a constant vector determined by initial conditions, and 

^ = -A- l T> (5.8) 
First, let us consider the subcase = 0. 

The solution of the system (5.3) (in the coordinate space), is provided by 
the following invariants: 
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A = . 1 {(At[c 2 + (b 2 + e 2 )cosuT])x + 
y/H + 1 w 



' ce r -, i b ■ \ 

— [—1 + cos lot\ suitor )y + 

"to 2 UJ 

be 6 

' [— 1 + cos lot] smujr)z}, (5.9) 

to to 



1 ce b 1 

A« = , {( — [-1 + coscjr] H smtor)x + ( — [c 2 + 

^/2a| + 1 ^ w cj 2 



(5 2 + e 2 ) costor])y + ( — [—1 + coscjt] smtor)z], (5.10) 

to 2 to 



1 r / &C r i e 

z/ = . ={( — 1 + coscjtJ H smu;r)x + 

v /2aITl - 



to* to 



^ r , n C 

— [— 1 + costurj H — sinwrjy + 



1 



4[6 2 + (c 2 + e 2 )coscjr])4, (5.11) 
'to 2 

where to = Vb 2 + c 2 + e 2 . 

From general formulas (5.9)— (5.11) we can obtain the following invariants, 
in a simple form: 

x = I^Td (cx - ev + bz) (5 ' 12) 



At = J ^J-^(ex + cy) 2 + b 2 (x 2 + y 2 ) + (c 2 + e 2 )z 2 + b(-2cxz + 2eyz) (5.13) 



_ (ex + cy) 2 + 6 2 (x 2 + y 2 ) + (c 2 + e 2 )z 2 + b(-2cxz + 2eyz) 

(cx — ey + bz) 2 
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The case where g a ^ 0, h Q ^ 0, r Q ^ can be obtained from the previous 
one by making the change of variables: 

a 2 g + c(cg - e/i p + br ) - a(6/i p + er Q ) 
x — > rr H — — , ^o.loj 

a{a z + b z + c z + e z ) 

t abg D - ceg Q + (a 2 + e 2 )h a - (ac + 5e)r Q 
a(a z + b z + c z + e z ) 

bcg Q + aeg + (ac - be)h + (a 2 + 6 2 )r p 
a(a z + b z + c z + e z ) 



Case II) a = 0. 

For r* = we obtain the following invariants: 

A = (cx-ey + bz) (5.18) 
fi = x 2 + y 2 + z 2 (5.19) 

V6 2 + c 2 + e 2 ^ - tt arctan 63 \ + ^ — (5.20) 

a b{—cx + ey) + {c z + e A )z 

Notice that A represents an invariant plane, whereas fi is an invariant 
spheric surphace centered at the origin. 
For g Q 7^ 0, h a ^ 0, r Q ^ 0, we obtain: 



A = cx - ey + bz - ^(c#o - eh + br ) 



(5.21) 



1 t be be 

cos(w-;( — + -»/io + 



1 t ,e 

smu-A—go 



Vc^Te^ 



d /v cu 2 ' 
dV 



c 2 + e 2 
~J 2 ~ 



r — ex — ey) + 



c , be be c 2 + e 2 . 

-h x + —y + z) (5.22) 

UJ UJ U LO 
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e c bcg - beh - x r 2 be be v . eg + ch 2 
z/ = (-z + -y 5 ) + ( x H y H — 2 H ) . 

X X v X ujx uj lux 

(5.23) 



Eq. (5.23) represents a cilinder with radius 

r = |A| 

and axis parallel to vector 

-x y + -z (5.24) 

u; id id 

passing throw the point of coordinates: 

2b c e go + b (c 2 — e 2 ) h — e x 2 m gob (c 2 — e 2 ) — 2bceh — cr x 2 



99 ' 99 ' 

- ^ + C/l ° (5.25) 

In the space of velocities, the invariant curves, which are now functions 
of ui,u 2 ,u 3 and p,can be deduced by the previous ones by performing the 
substitution 

a — > —a, x — > Mi, y — > 1^2, z — > 1*3, t — > p, d — > /c . (5.26) 

At this point let us examine some particular cases which emerge by de- 
manding that some arbitrary constants appearing in (5.2) are vanishing. Pre- 
cisely, we shall distinguish the cases: 

Subcase I: £ 4 7^ 0, where 



II 
12 
13 
14 



d ^ 0, a = b = c = e = 0, 

d^O, a = b = e = 0, c^O, 

d ^ 0, a = b = 0, c ^ 0, e ^ 0, 

d ^ 0, a = 0, b ^ 0, c ^ 0, e ^ 0. 



Subcase II: £ 4 = 0, with a — d — 0, and 
III : 6 = e = c = 0, 
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112 : b ^ 0, e = c = 0, 
113: b^O, c^O, e = 0, 
114 : b ± 0, c ^ 0, e ^ 0. 

The expressions for the invariants A, x an d ip corresponding to the sub- 
cases II,. ..,14 and III,.. .,114 are given by 

Subcase II 



Subcase 12 



Subcase 13 



\ = t- —x, (5.27) 
9o 

X = V--x, (5.28) 
9o 

t/; = 2--i. (5.29) 
9o 



\ = t- —x, (5.30) 
9o 

h + cz cx cy-r Q . cx 
X = cos 1 sin — , (5.31) 

c 9o c g Q 

cy-r cx h + cz . cx 
ip = cos sin — . (5.32) 

c 9o c g 



\ = --x + y, (5.33) 

e 

X = x 2 + y 2 + z 2 , (5.34) 



, d z\/ c 2 + e 2 . . 

ip — t-\ arctan . (5.35) 

Vc 2 + e 2 ex + cy 
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Subcase 14 



cx — ey 

X = r^ + z, (5.36) 



X = x 2 + y 2 + z 2 , (5.37) 

d b 2 x + eiex + cy) — bcz 

ib = t arctan -== v y> . 5.38 

^2 + c 2 + e 2 y/p + f + e^by + ez) v ' 



Subcase III 



Subcase 112 



Subcase 113 



X = t, (5.39) 



x = y--x, (5.40) 
9o 

^ = z - r ^ x . (5.41) 
9o 



X = t, (5.42) 



y = (x — ^) cos —z — (y + ^r) sin —2, (5.43) 

v 6 ' r D v 6 7 r Q v 7 

/ h Q . . b g . b 

tp — (x — — ) sin — z + (y + — J cos — z. (5.44) 
b r Q b r Q 



\ = t, (5.45) 



X = z + ^x, (5.46) 
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ijj = x 2 + y 2 + z 2 . 



(5.47) 



Subcase 114 



X = t, 



(5.48) 



X = cx — ey + bz, 



(5.49) 



; 2,2,2 
1p — X + y + Z . 



(5.50) 



6 Conclusions 



We summarize the main results achieved in this work. Using a group- 
theoretical framework, we have established a connection between the sym- 
metry generators of the incompressible 3D Navier-Stokes equations and the 
existence of solutions having a vortical dynamics. Precisely, by means of a 
procedure of symmetry reduction of Eqs. (1.1)— (1.2), we have found some 
novel (at the best of our knowledge) exact solutions which contain as special 
cases important solutions of the vortex-type well-known in the literature, 
as the Burgers vortex and shear-layer solutions. Such solutions have been 
discussed by the point of view of their possible physical meaning. 

It is suitable to recall that an approach of the algebraic and group- 
theoretical type to the equations for the fluid dynamics, especially to the 
incompressible Navier-Stokes equations, is not new and has been carried 
out by many authors (see, for example, [10], and [17]-[22]). In particular, 
the case of two space variables have been investigate more frequently and 
in more detailed manner. This is due to the fact that in 2D dimensions 
the possibility of defining a stream function [17] simplifies the mathematical 
treatment of the problem. In this spirit, recently Ludlow, Clarkson and Bas- 
som [17] have obtained non-classical symmetries for the 2D Navier-Stokes 
equations applying the so-called direct method [23]. On the other hand, 
in [20] Fushchich, Shtelen and Slavutsky have determined solutions of the 
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3D Navier-Stokes equations performing a systematic study of the inequiva- 
lent ansaetze of codimension 1 by means of which a direct reduction of the 
equations under consideration to ordinary differential equations is obtained. 

In our paper, together with the finding of exact solutions linked to the 
reduction method and together with the analysis of the role that the symme- 
try infinitesimal operators plays in understanding the vortical phenomena, 
we have classified in the context of the group theory all the invariant surfaces 
related to boundary conditions of physical interest. This kind of research is 
of special utility in possible applications to problems where the temperature 
is present. In particular, we mention the case in which to the Eqs. (1.1)— (1.2) 
one adds a transport equation describing the behavior of the temperature. 
These equations are known as the Navier-Stokes-Fourier equations, where 
the real coupling among them occurs just via the boundary conditions [24]. 

To conclude, we notice that in the symmetry algebra commutators as- 
sociated with the 3D Navier-Stokes Eqs. (1.1)— (1.2) arbitrary functions of 
integration appear. This indicates that the algebra is infinite-dimensional. 
Hence, it should been interesting to explore deeply this algebra and to clas- 
sify, for example, its one-dimensional subalgebras in order to discover other 
possible solutions of the 3D Navier-Stokes equations. 

Appendix A: the symmetry algebra 

The appearance of arbitrary functions in the vector fields (2.7)-(2.15) en- 
sures that the symmetry algebra of the Navier-Stokes equations is infinite- 
dimensional. This algebra is defined by the commutation relations 

[Vi(gi), V 1 (g 2 )] = V 4 (g 2 g\ - 9l g 2 ), [V 2 (h 1 ),V 2 (h 2 )] = V 4 (h 2 h x -h x h 2 ), 

(A.l) 

[Vsin), V 3 (r 2 )} = V 4 (r 2 f\ -n r 2 ), [V 4 (h), V 4 (k 2 )\ = 0, (A.2) 

[V^g), V 2 (h)\ = 0, [Vi(g), V 3 (r)\ = 0, (A.3) 
[Vi(g), V 4 (k)} = 0, [V^g), V 5 (a)} = V 1 (ag - 2a 9 t), (A.4) 
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[Vi(g), V 6 (b)\ = -V 2 (bg), [V 1 (g), V 7 (c)\ = 0, (A.5) 

[V^g), V 8 (d)\ = -V 3 (dg), [V^g), V 9 (e)\ = -V^e g), (A.6) 

[V 2 (h), V 3 (r)\ = 0, [V 2 (h), V 4 (k)\ = 0, (A.7) 

[V 2 (h), V 5 (a)\ = V 2 (ah -2aht), [V 2 (h), V 6 (b)] = V^bh), (A.8) 

V 2 (h), V 7 (c)} = -V 3 (ch), [V 2 (h), V 8 (d)} = 0, (A.9) 

[V 2 (h), V 9 (e)} = -V 2 (e h), [V 3 (r), V 4 (k)} = 0, (A.10) 

[V 3 (r), V 5 (a)} = V 3 (ar - 2a r t), [V 3 (r), V 6 (b)} = 0, (A.ll) 

[V 3 (r), V 7 (c)} = V 2 (cr), [V 3 (r), V 8 (d)} = V 1 (dr), (A.12) 

[V 3 (r), V 9 (e)} = -V 3 (e r), [V 4 (k), V 5 (a)\ = -2V 4 (ak + akt), (A.13) 

[V 4 (k), V 6 (b)} = 0, [V 4 (k), V 7 (c)} = 0, (A.14) 

[V 4 (k), V 8 (d)} = 0, [V 4 (k), V 9 (e)} = -V 4 (e k), (A.15) 

[V 5 (a), V 6 (b)] = 0, [V 5 (a), V 7 (c)} = 0, (A.16) 

[V 5 (a), V 8 (d)} = 0, [V 5 (a), V 9 (e)} = -2V 9 (ae), (A.17) 

[V 6 (b), V 7 {c)\ = -V s (bc), [V 6 (b), V 8 (d)} = V 7 (bd), (A.18) 

[V 6 (b), V 9 (e)} = 0, [V 7 (c), V 8 (d)} = -V 6 (cd), (A.19) 

[V 7 (c), V 9 (e)} = 0, [V 8 (d), V 9 {e)\ = 0. (A.20) 

Scrutinizing these commutation rules, it emerges that they contain the 
subalgebras 

{Vi,v 2 , y}, {v x y 2 , v 6 }, {v ± y 2 , v 9 }, {v ± y, y}, {v lt v 3 , v 8 }, {v lt v 3 , v 9 }, 
Wy 4 , y}, {v^y, v 9 }, {v 2 y, y}, {v 2 y, y}, {v 2 y 3 , v 9 }, {v 2 y 4 , y}, 
{v 2 y 4 , v 9 }, {v 2 y, v 9 }, {v 3 y 4 , y}, {v 3 y 4 , v 9 }, {yy, v 9 }, {v 4 y, v 9 }, 

{yy,y}, 

(A.21) 

which turn out to be isomorphic to the algebra of the Euclidean group E 2 . 
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Appendix B: derivation of the solution (4.41a)— 
(4.41d) 



In order to find a particular solution of Eqs. (4.41a)-(4.41d), let us assume 
that 

A 1 = h(6),A 2 = f(6),A 3 = g(9). (B.l) 
Then, Eq. (4.41d) provides 

Ai = c0*, (B.2) 
c being a constant of integration, while Eqs. (4.41a)-(4.41c) give 

1 3 

2tt - 29n e - r]Ti v - £tt 5 + -c9^ = 0, (B.3) 

H = ^F(0), (B.4) 

7r„ = - 2 G{6), (B.5) 

respectively, where 

F(6) = -{2v + c6^)f + (2i/<9 + 2c#i - 6 2 )f - 4v6 2 f", (B.6) 

G{0) = -{2u + + (2v6 + 2c^i - 6 2 )g' - 4v6 2 g", (B.7) 

with f =%,</=%. 

By integrating Eq. (B.4) with respect to £, we get 

7r = -^F(0) + a(77,0), (B.8) 

where 01(77, Q) is a function of integration. 
From (B.8) and (B.5) we have 

TTr, = a v = (B.9) 
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from which 

a = - C ^l + p(0), (B.10) 

where (3 is an arbitrary function of 9. Taking account of (B.10), Eq. (B.8) 
becomes 

K= _m_m +m (B11) 

Then, substitution from (B.ll) into Eq. (B.3) yields 

F = a9'i, (B.12) 



G = b6*, (B.13) 

P=^ + co0, (B.14) 

where a, b, Co are constants. 

Now it is convenient to put 9? = p. So, with the help of (B.12), Eq. 
(B.6) can be written as 

vp 2 f PP + (\p 3 - cp 2 - 2up)f p + (cp + 2u)f + ap 3 = 0. (B.15) 
The general solution of this equation is 



f(p) = -lap + v^exp P(4 ^ P \ Cl + c 2 erf ^ J 20 ^* ]. (B.16) 

The form of g(p) can be determined in the same way. Coming back 
to the original variable 9 and exploiting (B.l), (B.2), (B.ll), the solution 
(4.41a)-(4.41d) is readily found. 
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